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Using first-principles calculations, the electronic and magnetic properties of orthorhombic BaFeOs 
(BFO) are investigated with local spin density approximation (LSDA). The calculations reveal that 
at the optimized lattice volume BFO has a lower energy in ferromagnetic state as compared with 
antiferromagnetic state. At the equilibrium volume, BFO shows metallic behavior, however, under a 
large tensile strain (~ 25%), BFO shows half-metallic behavior consistent with the integer magnetic 
moment of 4.0/iB/fu mainly caused by the t 2 g and Cg electrons of Fe. Including a Hubbard-like 
contribution U (LSDA+[/) on Fe d states induced half-metallic bahvior without external strain, 
which indicates that U can be used to tune the electronic structure of BFO. The magnetic moments 
remained robust against ~ 10% compressive and tensile strain. At large compressive (tensile) strain, 
the half-metallicity of BFO is mainly destroyed by the Fe-d (O-p) electrons in agreement with the 
non-integer value of the magnetic moments of BFO. 


PACS numbers: 71.15.Mb,71.27.-|-a,71.30.-|-h,77.80.bn 


Since the discovery of half-metallic materials in 1983 
by Groot et.al^ many other new materials have been dis¬ 
covered using first-principles calculations because these 
materials are promising candidates in the area of spin- 
electronicsj^ Perovskite-structure oxides, ABO3, are of 
technological interest due to their numerous and of¬ 
ten coupled functionalities, such as ferroelectricity, mag¬ 
netism, multiferroicity, metal-insulator transitions, etc. 
The perovskite family is known to have a large vari¬ 
ety of properties'^— including spin dependent transport 
properties^. Some perovskite materials also show half- 
metallicity driven by external pressure. When the exter¬ 
nal pressure is below ^ 80GPa,NiGrO3 is an antiferro¬ 
magnetic (AFM) semiconductor and behaves as an AFM 
half-metal when the external pressure is increased^. At 
ambient conditions Sr2NiRe06 is a semiconductor, and 
through external pressure, a semiconductor to half-metal 
transition is observed— La2GoMn06 also transformes to 
half-metallic state under external pressurei^. 

Fe-containing ferroelectric oxides are also attracting 
tremendous attention due to their magnetic properties. 
Particularly, those oxides that exhibit magnetic and fer¬ 
roelectric characteristics simultaneously can have practi¬ 
cal device applications such as spin transistor memories, 
whose magnetic properties can be tuned by electric fields 
through the lattice strain effect— i Fe-based perovskite 
BaFe03 (BFO) normally assumes a hexagonal crystal 
structure, although various polymorphs have been ob¬ 
served with oxygen deficiencyi^^— Recent experiments 
observed, however, ferromagnetism in cubic BFO, where 
the metastable cubic structure was shown to be stable 
down to 8K— ^ Very recently, Li et al— ^ used DFT-bG to 
investigate the helical structure of cubic ferromagnetic 
(FM) BFO and found AFM to FM transitions under ex¬ 
ternal pressure. 

It has been shown that the functional behavior in ox¬ 
ides can be enabled or enhanced when they are prepared 
in the form of thin films. For example, ferroelectricity 
has been reported in thin hlms of SrTiOs, which is a 
quantum paraelectric in its bulk form, and enhanced po¬ 


larization has been observed in thin films of ferroelec¬ 
tric BaTiOs— In both cases, the improved behavior 
is attributed to the bi-axial epitaxial strain introduced 
when the film is grown coherently on a substrate with 
a different lattice constant. As there are no theoreti¬ 
cal and experimental work have been performed on or¬ 
thorhombic BaFeOs, and the possibility of orthorhombic 
phase can not be ignored when we synthesis BaFeOs. 
Such orthorhombic phase is already exists in other per- 
ovskites, e.g., BiFeOs, BaTiOs, etc. Therefore, we fo¬ 
cus on orthorhombic BaFeOs and show that orthorhom¬ 
bic BaFeOs has a negative formation energy, which indi¬ 
cates that we can not ignore this crystallographic phase 
of BaFeOs. In particular, we also show that that exter¬ 
nal strain and strong correlation can induce half-metallic 
ferromagnetism in orthorhombic BaFeOs. 

The First-principles calculations based on density 
functional theory are performed with plane-wave and 
pseudopotential method as implemented in Quantum 
Espresso package—^ The exchange correlation effects are 
treated within local spin density approximation(LSDA). 
Later on, the on-site Goulomb potential [/(=5.0eV)^ 
has also been added in LSDA to perform LSDA-I-C/ 
calculations to correctly describe the electronic struc¬ 
ture of orthorhombic BaFeOs (denoted as BFO in the 
manuscript). The ultrasoft pseudopotentials are used to 
describe the core-valence interactions. The valence wave 
functions and the electron density are described by plane- 
wave basis sets with kinetic energy cutoff of 30 Ry with 
12 X 12 X 12 Monckhorst-Pack grid. All the computa¬ 
tional parameters were carefully converged, and the con¬ 
vergence criterion for the electronic energy is 10“® eV. 

It is essential to optimize the lattice volume of BFO 
either using the molecular dynamics or the total energy 
calculations. We followed the later procedure and calcu¬ 
lated the total energies in the nonmagnetic (NM) and fer¬ 
romagnetic (FM) states. Using the equation of state^i, 
the equilibrium lattice volume Vq is determined, which 
is 56.28 A^/formula unit (fu) (a = 3.82,6 = 5.43, c = 
5.45A). Figure UKa) shows the strain, which is , vs 
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FIG. 1: (a) The LSDA calculated total energy/fu vs strain 
for FM(filled circles) and NM(empty circles) BFO. (b) The 
calculated total magneitc moment/fu in LSDA(filled circles), 
LSDA+[/(empty circles), and the local magnetic moment of 
Fe in LSDA(triangles) vs strain of BFO. The magnetic mo¬ 
ments of BFO are shown in the units of ^ b - The crystal struc¬ 
ture of orthorhombic BFO is also shown, where cyan, red, and 
green balls represent Ba, Fe, and O atoms, respectively. 

energy, where positive (negative) strain corresponds to a 
tensile (compressive) strain. In the whole range of strain, 
the FM state remains in the lowest energy state. Us¬ 
ing the optimized lattice volume, the formation energy 
of BFO is calculated^^, which is ^ —3.80 eV/fu. The 
negative formation energy illustrates that it may not be 
difficult to synthesize BFO. To see any possible antiferro¬ 
magnetic (AFM) interactions between the Fe atoms, the 
FM and AFM interactions between the Fe atoms were 
also considered. We found that the FM state is 0.130 
eV/fu lowered than the AFM state. It should be men¬ 
tioned that the non-stoichiometric BFO may have other 
magnetic structures.— Once we have confirmed that the 
FM state is the ground state of BFO, then we inves¬ 
tigated the magnetic and electronic properties of BFO 
in the FM state under different tensile and compressive 
strains. 

The calculated magnetic moment/fu is shown in 
Fig. [TJb). It is evident that the magnetic moment of 
the lattice is very sensitive to the strain, and the mag¬ 
netic moment decreases with compression. The magnetic 
moment increases monotonically with the tensile strain 
up to 25%. The increment/decrement of the mag¬ 
netic moment is mainly due to elongation/compression of 
the Fe-0 bond length which further decreases/increases 
the hybridization between the Fe-d and 0-p orbitals that 
control the charge transfer between the Fe and O atoms. 
Beyond ^ 25% tensile strain, saturation in the magnetic 
moments, 4.0/rB /fu, i.e., Fe^+ can be seen. The local 
magnetic moments of Ba, Fe, and O were also analysed, 
and the increase/decrease in the total magnetic moments 
is mainly contributed by the local magnetic moments of 
the Fe atoms, which are also shown in the same figure. 


Small dependence of the O magnetic moments on strain 
is also observed. The LSDA calculated total magnetic 
moment of BFO at the equilibrium volume is 3.34/rB/fu, 
whereas the local magnetic moments of Fe and O are 
2.74/xb and 0.19/rB- The integer (saturated) magnetic 
moment at a large tensile strain is an indication of half- 
metallic ferromagnetism in BFO. Therefore, we can ex¬ 
pect half-metallic ferromagnetism in the strained BFO, 
which will be discussed in the following paragraphs. The 
local magnetic moment of Fe (O) at ^ 25% tensile strain 
is 3.35/rB(0.20/rB) which further shows that the Fe local 
magnetic moment is mainly responsible for the integer 
magnetic moment of BFO. The magnetic interactions be¬ 
tween the Fe and O atoms remained ferromagnetic due 
to positive induced polarization at the O sites. 

We observed that the ^ 25% strained BFO can show 
interesting magnetic properties, and such a large strain 
may be achieved experimentally by depositing BFO on a 
large mismatched substrate. However, LSDA usually ig¬ 
nores the strong correlation effect between the Fe atoms. 
If such a strong correlation effect exists between the Fe 
atoms, then it is expected that BFO may show integer 
magnetic moment (half-metallicity) without strain. We 
then repeated the above calculations using the LSDA-I-U 
approach, and the calculated magnetic moment at each 
strain is also shown in Fig. [TJb). The total magnetic 
moment of BFO at Vo is 4.0/iB/fu, whereas the local 
magnetic moments of Fe and O are 3.73/rB and 0.09/iB- 
One can easily notice that including U has increased (de¬ 
creased) the magnetic moment of Fe (O). It is also in¬ 
teresting to note that the total magnetic moments are 
remained robust against the strain, and the magnetic 
moments are saturated to 4.0^B/fu in the whole strain 
region. The integer value of the magnetic moment is 
an indication of half-metallic magnetism in BFO, as we 
observed in the strained BFO LSDA calculations. How¬ 
ever, a drastic change in the magnetic moment can be 
seen at a large compressive or tensile strain, where the 
magnetic moment deviates from the integer value sug¬ 
gesting a possible non half-metallic behavior. Note that 
for ~ 25% tensile strained, LSDA showed half-metallic 
behavior and LSDA-I-C/ has destroyed the half-metallicity 
in BFO as shown below. Hence, there is a strong com¬ 
petition between the strain energy, which localizes the 
minority electrons (in case of elongation), and Coulomb 
energy that delocalizes electrons at large strain(in case of 
compression), and strain energy looses at the cost of half- 
metallicity. From these observations, it can be inferred 
that a large strain can destroy the half-metallic nature of 
BFO even if the effect of U is included, which will also 
be highlighted in the following paragraphs. 

As shown above that BFO can have integer mag¬ 
netic moment (expected half-metalic ferromagnetism) at 
a large strain in the LSDA calculations, however, BFO 
shows integer magnetic moment at equilibrium volume 
and a large strain can suppress the integer magnetic mo¬ 
ment (expected non half-metalic ferromagnetism) nature 
of BFO in the LSDA-I-U calculations. To further confirm 
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FIG. 2: (Color online) The calculated total electronic den¬ 
sity of states of BFO under different strain in LSDA(a) and 
LSDA-l-17 (b). All the studied strain are mentioned where 
0.0% strain corresponds to the equilibrium structure of BFO. 
The Fermi energy Ep is set to zero. The upper (lower) panels 
shows majority (minority)-spin DOS. 


these gestures, we calculated the total electronic den¬ 
sity of states (DOS) at each strain, and the total DOSs 
are shown in Fig. [5] (a), where 0.0% strain corresponds 
to the equilibrium volume. At the equilibrium volume, 
BFO is a metallic magnetic material with no trace of 
half-metallicity. The Fermi energy lies in the conduction 
band, and one can also see a pseudogap in the minority 
spins just below the Fermi energy. The pseudogap can 
be shifted to the Fermi energy, if an external strain is 
imposedi^ By a compressive strain, BFO shows metallic 
behavior and the pseudogap is shifted more below the 
Fermi energy, i.e, the Fermi energy penetrates more the 
conduction band. Therefore, the compressive strain does 
not help in inducing half-metallicity. However, on the 
other hand, for a tensile strain the Fermi energy shifts 
towards the pseudogap, and the Fermi energy shifts away 
from the conduction band, and at a large strain ( 25%) 
the Fermi level lies just at the conduction band edge Ec 
indicating that the minority spin band is completely oc¬ 
cupied. Such movement of the Fermi energy in the mi¬ 
nority spins state increases the magnetic moment [see 
Fig. mb)]. Beyond 25% tensile strain, the Fermi energy 
clearly separates the conduction and valance bands and 
a bandgap of ~ 0.07 eV in the minority spins state can 
be seen. Such an electronic structure of a FM material 
is known as a half-metal ferromagnet consistent with the 
integer magnetic moment as discussed above. 

We have shown that strain can easily tune the elec¬ 
tronic structure of BFO. However, for oxides contain¬ 
ing transition metals, the strong on-site Coulomb energy 
and exchange interactions for the d electrons between 
neighbouring transition metal ions due to the local d 
orbitals can not be ignored. Knowing that LSDA usu¬ 
ally underestimates the size of the bandgap in a mate¬ 
rial with strongly localized d electrons and even predicts 
the metallic behvior of meterials that are known to be 
insulator. We then used LSDA-I-U to describe the cor¬ 
rect electronic structures of BFO under strains, which 
are shown in Fig. [2l[b)- The band structure at the opti¬ 
mized volume is represented by 0.0% strain, and one can 
easily see a large bandgap 1.04 eV) in the minority 
spins state, and the Fermi level lies in the bandgap-a 



FIG. 3: (Golor online)The LSDA calculated total(lower panel) 
and partial density of states(upper panel) of BFO under - 
16.50% (a), equilibrium (b), and 25.00% (c) strain. In the up¬ 
per panel, the solid(black), dashed(blue), and dotted(green) 
lines shows Fe-t 2 g, Fe-Cg, and O-p states, respectively. The 
Fermi energy Ep is set to zero. The upper (lower) panels 
shows majority (minority)-spin DOS. 


half-metal ferromagnet. These calculations indicate that 
the strong Coulomb interaction between the d electrons 
has overcome the tensile ~ 25% strain that was required 
for half-metallic ferromagnetism in BFO in the LSDA 
calculations. The main effect of strain is to shift the 
Fermi energy towards valence band edgeEy in the mi¬ 
nority spins state in the strained BFO, and for a large 
tensile strain Ep crosses Ey. Thus some of the minor¬ 
ity spins state becomes unoccupied due to a large ten¬ 
sile strain that further increases the magnetic moment of 
BFO and no longer remains integer magnetic moment. 
And the compressive strain shifts the Fermi energy to Ey 
that promotes more electrons from the conduction band 
to the valance band in the minority spins state and the 
magnetic moment decreases from its integral value. Thus 
a large compressive or tensile strain can destroy the half- 
metallic nature of BFO when the effect of U is considered. 
Therefore, the main effect of U is to push apart the bond¬ 
ing and anti-bonding states that widen the bandgap, but 
strain shifts the Fermi energy and can drastically change 
the electronic structure. Figure [ll[b) also demonstrates 
that BFO can retain its half-metallic nature under ± 10% 
strain, and beyond this strain BFO will have a metallic 
nature. The half-metallic ferromagnetic behvior is con¬ 
sistent with the integer magnetic moment as discussed in 
Fig. mb) and the decrease/increase in the magnetic mo¬ 
ment at a large compressive/tensile strain is due to the 
non half-metallic property of BFO. 

It is very essential to further investigate the atomic 
origin of strain or correlation induced half-metallic ferro¬ 
magnetism in BFO, and why BFO looses its half-metallic 
nature in the presence of a large strain even when U is 
included. Figure [3] shows the LSDA total and atom pro¬ 
jected orbital resolved DOS of —16.5% (panel (a)), 0.0% 
(panel (b)), and 25% (panel (c)) strains. For 0.0% strain, 
the partial DOS shows strong p-d bonding which gives 
rise to metallicity in BFO. The Fermi energy crosses the 
conduction and valence bands, and the minority spins 
band near the Fermi energy is mainly dominated by 
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FIG. 4: (Color online)The LSDA+17 calculated total(lower 
panel) and partial density of states(upper panel) of BFO 
under -16.50%(a), equilibrmm(b), and 16.85%(c) strain.In 
the upper panel, the solid(black), dashed(blue), and dot- 
ted(green) lines shows Fe-t 2 g, Fe-eg, and 0-p states, respec¬ 
tively. The Fermi energy Ef is set to zero. The upper (lower) 
panels shows majority (minority)-spin DOS. 


the Fe eg and t 2 g electrons. The degenerate behavior 
of Fe eg and t 2 g indicates that the crystal field is not 
strong enough to remove the degeneracy of the d or¬ 
bitals. The tensile strain mainly narrows the d band, 
and the bandwidth of the d band is also decreased. Due 
to strong overlap of the Fe-d and 0-p-orbitals, the O- 
PDOS also becomes localized and participates in the half- 
metallicity of BFO. For the compressive strain, the bands 
become very wide, and these are basically the minority 
t 2 g and eg bands that help in producing/destroying the 
half-metallicity in BFO. The LSDA-bD calculated PDOS 
of—16.5%(panel (b)), 0.0% (panel (b)), and 16.85%(panel 
(c)) strains are also analysed in Fig. 0] One can clearly 
see that at the equilibrium volume, the half-metallicity 
is mainly caused by the shift of Fe t 2 g and eg electrons, 
and these bands are completely occupied in the majority 
spins state, which were partially occupied without con¬ 
sidering U. In the minority spins state, the t 2 g and eg 
electrons of Fe are almost unoccupied, and one can see 
a large bandgap of 1.04eV due to the unoccupation of 
these bands. As discussed above that for the compres¬ 
sive strains, the Fermi energy goes deeper into the con¬ 
duction band, and now it is clear from Fig. SKa) that the 
half-metallicity is mainly destroyed by t 2 g and eg bands, 
and predominantly by the t 2 g electrons. However, for 
the tensile strain, the half-metallicity is not destroyed by 
the t 2 g and eg electrons of Fe, but by the p electrons of 
O. The tensile strain shifts the Fermi energy towards Ec 


and the top of the valance band in the minority spins is 
dominated by the p-orbitals of O and the Fe d band is 
far from the Fermi energy due to the on-site Coulomb 
repulsion(see Fig. SKc)). Finally, comparing the PDOS 
of LSDA with LSDA+[/, it is noticed that the hybridiza¬ 
tion between Fe-d and 0-p electrons is weaker in the 
LSDA-|-[/ case, and the top of the valance band in both 
the spins state is mainly dominated by the p electrons of 
oxygen. Before summarizing our theoretical calculations, 
we would like to comment on the helical magnetic order, 
which is observed in cubic BFOJ^ Magnetism is con¬ 
sidered to be a cooperative phenomenon and we mainly 
focused on the collinear behaviour of spin magnetic mo¬ 
ments of Fe in orthorhombic BFO. Investigating the he¬ 
lical magnetic order is beyond the scope of the present 
work. 

To summarize, we predicted strain and correlation in¬ 
duced half-metallic ferromagnetism in BFO using DFT 
calculations. A large strain would be needed for half- 
metallicity, but no strain is required to have half-metallic 
ferromagnetism, when the correlation effect between the 
Fe d electrons is considered. The half-metallic behavior 
in the strained BFO is mainly caused by the shift of the 
t 2 g and eg electrons of Fe from the Fermi energy. Large 
compressive and tensile strain in the presence of U can 
destroy the half-metallicity in BFO. The presence (ab¬ 
sence) of half-metallic ferromagnetism in BFO is consis¬ 
tent with integer (non-integer) magnetic moment of BFO. 
Therefore, from these DFT calculations, we can infer that 
U is playing a vital role in inducing half-metallic ferro¬ 
magnetism in BFO. Generally, it may be difficult to im¬ 
pose a large strain due to elastic modulus of a material 
and such a large strain may induce dynamic instability. 
Therefore, for materials with inherent correlation effect 
between the transition metals, e.g. Fe, can show interest¬ 
ing properties (half-metallic magnetism) with no external 
strain. Strain in a strongly correlated system can be used 
as a parameter to see the maximum absorption of the 
strain to retain its half-metallic behavior, and we found 
that BFO can withstand ~ ±10% strain before convert¬ 
ing into metallic BFO. As U can be considered as a source 
to overcome the large tensile strain, therefore, we believe 
that those materials that can have strong Coulomb inter¬ 
actions are good candidates for applications in the area 
of spin-electronics or multiferroics. Further experimental 
work is needed to confirm our predictions. 
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